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U(1)a Goldberger-Treiman relation is extended to accommodate SU(3) flavour breaking 
and the implications of a more precise measurement of the rj and ry'-nucleon couplings 
are discussed. A comparison with the existing literature on pseudoscalar meson decay 
constants using large- N c chiral Lagrangians is also made. 



"This research is supported in part by PPARC grants PP/G/O/2002/00470 and PP/D507407/1. 



1. Introduction 



The phenomenology of the pseudoscalar mesons opens a window on many interesting as- 
pects of the non-perturbative dynamics of QCD, including spontaneous chiral symmetry 
breaking, the electromagnetic and gluonic axial anomalies, the OZI rule, the gluon topolog- 
ical susceptibility, and so on. Nevertheless, until comparatively recently, phenomenological 
analyses did not take fully into account the role of the gluonic U (I) a anomaly in this sector, 
with the result that most of the existing determinations of the pseudoscalar meson decay 
constants and couplings are based on over-simplified theoretical formulae which miss the 
most interesting anomaly-sensitive physics. 

In a previous paper [1] (see also [2]) we have analysed in detail the radiative decays of 
the neutral pseudoscalar mesons, n ,rj,rf — ► 77, and provided a set of formulae describing 
these processes together with modified Dashen-Gell-Mann-Oakes-Renner (DGMOR) [3, 4] 
relations which fully include the effect of the anomaly and the gluon topological suscepti- 
bility in the flavour singlet sector as well as explicit flavour SU(3) breaking. In this paper, 
we confront these formulae with experimental data to extract a set of results for the four 
pseudoscalar meson decay constants, f 0r, ',f° v ,f 8r1 ' and f Sv . 

Our approach, which is based on a straightforward generalisation of conventional 
PC AC to include the anomaly in a renormalisation group consistent way, is unique in 
that it is theoretically consistent yet does rely on using the 1/N C expansion. The decay 
formulae and DGMOR relations we derive are valid for all N c . At some point, however, 
the fact that the 7/ is not a Nambu-Goldstone (NG) boson means that to make predictions 
from our formulae we need to augment the standard dynamical approximations of PCAC 
or chiral perturbation theory with a further dynamical input. Here, we make a judicious 
use of 1/N C ideas to justify the use of the lattice calculation of the topological suscep- 
tibility in pure Yang-Mills theory as an input into our flavour singlet DGMOR relation. 
The self-consistency of this approach allows us to test the validity of l/N c methods in the 
U(1)a sector against real experimental data. In particular, the flavour singlet DGMOR 
relation is the finite- -/V c generalisation of the well-known Witten-Veneziano formula [5, 6], 
and we can use our results to study how well the large- N c limit is realised in real QCD. This 
is especially important in view of the extensive use of large N c in modern duality-based 
approaches to non-perturbative gauge theories such as the AdS-CFT correspondence. 

The relation of our analysis of the low-energy decays 7/(77) — > 77 to high-energy physics 
in the form of the first moment sum rule for the polarised photon structure function gj [7, 8] 
is also discussed. The dependence of this sum rule on the photon momentum encodes many 
aspects of non-perturbative QCD physics and its full measurement is now within the reach 
of planned high-luminosity e + e _ colliders [9]. Our methods can also be readily extended 
to a variety of other reactions involving the pseudoscalar mesons, including rf (77) — > V7, 
where V is one of the flavour singlet vector mesons p,lo,4>, 7/(77) ~~ * 7r + 7r~7, a variety of 
electro and photoproduction reactions such as jp — ► pr]'(rj) and jp — ► p<p, and rj and 77' 
production in pp collisions, pp — > pprf(rj). 

The latter are relevant for determining the couplings g v NN and g v >NN which occur 
in generalisations of the Goldberger-Treiman relation [10]. In this paper, we present a 
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new version of the U{\)a Goldberger-Treiman relation [11, 12, 13] incorporating flavour 
octet-singlet mixing in the rj — rf sector. We are then able to use our results for the 
decay constants to show how the physical interpretation of this relation, which is closely 
related to the 'proton spin' problem in high-energy QCD (i.e. the first moment sum rule 
for the polarised proton structure function eft), depends critically on the experimental 
determination oi g v 'NN- 

Finally, in an appendix, we give a brief comparison of the similarities and differences 
between our approach and an alternative theoretically consistent framework for describing 
r\ and rf physics, viz. the large- N c chiral Lagrangian formalism of refs.[14, 15, 16] (see 
also [17, 18]). A comprehensive review of the phenomenology of the pseudoscalar mesons, 
incorporating the chiral Lagrangian results, is given in ref.[19] 

2. Radiative decay and DGMOR formulae 

In this section, we present the radiative decay and DGMOR formulae and discuss how to 
extract phenomenological quantities from them, emphasising their 1/N C dependence. A 
brief review of their derivation in the language of 1 U(1)a PCAC is given in section 3, but 
we refer to the original papers [1, 2, 20] for full details, including a more precise formulation 
using functional methods. Throughout, we consider the case of three light flavours, n/ = 3. 

The assumptions made in deriving these formulae are the standard ones of PCAC. 
The formulae are based on the zero-momentum chiral Ward identities. It is then assumed 
that the decay 'constants' (f aa {k 2 ) in the notation below) and couplings (^« 77 (k 2 )) are 
approximately constant functions of momentum in the range from zero, where the Ward 
identities are applied, to the appropriate physical particle mass. It is important that 
this approximation is applied only to pole-free quantities which have only an implicit 
dependence on the quark masses. This is equivalent to assuming pole-dominance of the 
propagators for the relevant operators by the (pseudo-)NG bosons, and naturally becomes 
exact in the chiral limit. Notice that these dynamical approximations are equally necessary 
in the chiral Lagrangian approach, where they are built in to the basic structure of the 
model - the effective fields are chosen to be in one-to-one correspondence with the NG 
bosons and pole-dominance is implemented through the low-momentum expansion with 
momentum-independent parameters corresponding to the decay constants. 

The relations for the 7r° decouple from those for the 77 and rf in the realistic approxi- 
mation of SU(2) flavour invariance of the quark condensates, i.e. (uu) = (dd). The decay 
formula for ir° — > 77 is the standard one, 



7T<?7T77 — ®< 




(2.1) 



with <2p, 



'cm 



3 N c , while the DGMOR formula is simply 



final = - (m u (uu) + m d {dd)) 



(2.2) 



The coupling g nll is defined as usual from the decay amplitude: 

( 77 |7T°) = - i g nll e Xpa pPiP2e X (pi)e p {p 2 ) 



(2.3) 
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in obvious notation [1]. The r.h.s. of eq.(2.1) arises from the electromagnetic U(1)a anomaly 
and is of course sensitive to the number of colours, N c . The decay constant has a simple 
theoretical interpretation as the coupling of the pion to the flavour triplet axial current, 
(0| J^ 5 |7r) = ikfj,f w , which is conserved in the chiral limit. 

In the 77, 7/ sector, however, explicit SU(3) flavour breaking means there is mixing. 
The decay constants therefore form a 2 x 2 matrix: 

/ fOt]' fOr] \ 

f aa = [fsn' jsn) (2-4) 

where the index a is an SU(3) x U(l) flavour index (including the singlet a = 0) and a 
denotes the physical particle states 7r°, 77,77'. These four decay constants are independent 
[14, 15, 21, 22], in contradiction to the original phenomenological parametrisations which 
erroneously expressed f aa as a diagonal decay constant matrix times an orthogonal 77 — 77' 
mixing matrix, giving a total of only three independent parameters. Sometimes the four 
decay constants are expressed in terms of two constants and two mixing angles but, while we 
also express our results in this form in section 4, there is no particular reason to parametrise 
in this way. 

The other new feature in the 77, rf sector is the presence of the gluonic U(1)a anomaly 
in the flavour singlet current. This means that even in the chiral limit the 77' is not a true 
NG boson and therefore the direct analogue of the PCAC formula (2.1) does not hold. 
Nonetheless, we can still write an analogous formula but with an extra term related to the 
gluon topological susceptibility. 

The decay formulae are 1 : 

/° T/ <V77 + /°% 77 + y/6A gGry = a° cm ^ (2.5) 

/ ^ Sjj'77 + / ^drfyy "em (2-6) 

where a° m = |^^c and af m = t^N c , and the corresponding DGMOR relations are: 

(/(V)2 m 2 ; + (/0^2 m 2 = _ 2( mu (uu) + m d {dd) + m s {ss)) + 6A (2.7) 

ftofftof m n + f°vfSv m 2 = _ ^(^(uu) + m d (dd) - 2m s (ss}) (2.8) 

(/ 8 "') 2 m; 2 + (/ 8 ") 2 m2 = - \(m u {uu) + m d {dd) + Am s {ss)) (2.9) 



1 Notice that compared to refs.[l, 2, 20, 9] we have changed the normalisation of the flavour singlets by 
a factor l/y^n/ = l/\/6. The normalisations here are txT a T b — \5 a b for generators T a (a = 0, 1, . . . 8) 
of 51/(3) x U(l), i.e. T l = |A 4 (i = 1, . . . 8) where A 1 are the Gell-Mann matrices, and T° = 4-1 for the 
flavour singlet. The d-symbols are defined by {T a ,T b } — d abc T c and include dooo = do33 = doss = 0(330 = 

dsSO = V/|; ^338 = ^383 = —^888 = v/|- 
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There are several features of these equations which need to be explained. First, because 
of the anomaly, the decay constants in the flavour singlet sector can not be identified as 
the couplings of the 77, rf to the singlet axial-vector current. In particular, the object F®, 
defined as (0\J® 5 \rf) = ik^F®, is not a renormalisation group invariant and so is not to be 
identified as a physical decay constant. The derivation of the radiative decay and DGMOR 
formulae in the next section makes it clear that F®, plays no role - certainly it is not 
f 0f i'. Instead, our decay constants are defined in terms of the couplings of the it, rj, rf to 
the pseudoscalar currents through the relations f aa (0\<p^\^ a ) = d a i) C ((f) c ) (for notation, see 
section 3). This coincides with the usual definition except in the flavour singlet case. 

The coefficient A appearing in the flavour singlet equations is the non-perturbative 
constant that determines the topological susceptibility in QCD. Recall that the topological 
susceptibility x(0) is defined as 

X(0) = J d 4 x i(0\T Q(x) Q(0)|0) (2.10) 

where Q = ^tiG^G^ is the gluon topological charge. The anomalous chiral Ward 
identities determine the dependence of x(0) on the quark masses and condensates up to a 
single non-perturbative parameter [23, 1], viz: 

or alternatively, 

—A m u m ( irn s {uu){dd){ss) 



X(0) 



m u nidm s (uu) (dd) (ss) — A(m u md(uu) (dd) + m u m s (uu)(ss) + m ( im s (dd)(ss)) 

(2.12) 



Notice how the well-known result that x(0) = if any of the quark masses is zero is realised 
in this expression. 

The final new element in the flavour singlet decay formula is the coupling parameter 
gc-y-y- This is unique to our approach. It takes account of the fact that, because of mixing 
with the pseudoscalar gluon operator Q due to the anomaly, the physical rf is not a NG 
boson. 5G77 is n °t a physical coupling, although it may reasonably be thought of as the 
coupling of the photons to the gluonic component of the rf . (A justification for this picture 
is given later.) Its precise field-theoretic definition is given in section 3 and refs.[l, 2, 20]. 
We should emphasise, however, that an analogous parameter must appear in any correct 
current algebra analysis involving the 7/ , such as the 1/N C chiral Lagrangian approach of 
refs.[14, 15, 16]. Phenomenologically, it may be found to be relatively small (we test this 
numerically in section 4), but any formulae which leave out such a term completely are 
inevitably theoretically inconsistent. 

Of course, in addition to these formulae, we also have the DGMOR relations for the 
remaining pseudoscalars. For example, for the K + , 

fK m K = - (m u {uu) + m s {ss)) ( 2 -13) 



- 4 - 



If we assume exact SU{2) flavour symmetry, we can then use eqs.(2.2) and (2.13) to elim- 
inate the quark masses and condensates in favour of f n ,fK,m 2 and m 2 K in the DGMOR 
formulae for the rj and rf. We find: 

(/°"') 2 m2, + (/0^ m 2 = l(/2 m 2 +2/ 2 rm 2 :)+6i4 (2 . 14) 

fXfVmZ + ftftoirn* = ^#(/>* - f K m\) (2.15) 

(/ V)2 m « + (/ a J) 2 m 2 = _l U 2 m 2_ Af 2 Km 2 K) (216) 

A quick look at the full set of decay and DGMOR formulae (2. 5), (2. 6) and (2.14)-(2.16) 
now shows that we have five equations for six parameters, viz. f° v ,f° v ,f 8ri ,f 8v ,A and 
gcyy, assuming that f n , fx and the physical masses are known along with the experimental 
values for the couplings 5y 77 and grfyj. It is not surprising that this set is under-determined. 
In particular, the necessary presence of the unphysical coupling gcyy in the flavour singlet 
decay equation essentially removes its predictivity. At this stage, the best we can do is 
to evaluate the singlet decay constants and the gluonic coupling gc-yy as functions of the 
topological susceptibility parameter A. This is done in section 5. In order to make more 
progress, we therefore need a further, dynamical input. 

As yet, everything we have done has been entirely independent of the 1/N C expansion. 
The 1/N C expansion (or the OZI limit - see ref.[24] for a careful discussion of the differences) 
is known to give a good approximation to many aspects of the dynamics of QCD and could 
provide the required extra input, but its application to the U(1)a sector needs to be handled 
with great care. For example, in the chiral limit, the mass of the rf, which arises due to the 
anomaly, is formally m 2 , = 0(1/N C ); however, numerically (allowing for the quark masses 
going to zero) this is of the same order of magnitude as a typical meson such as the p and is 
certainly not small. Generally, it is not clear that quantities which are formally suppressed 
in 1/N C are in fact numerically suppressed in real QCD, so we must be extremely careful 
in applying 1/N C methods here. 

Despite these caveats, we will see that 1/N C can play a useful role in analysing the 
decay formulae and DGMOR relations. Conventional large N c counting gives the following 
orders for the various quantities: f aa = 0(s/N c ) for all the decay constants; the rj and 
rf couplings g^j-y = 0(y/N c ) but gc-y-y = 0(1); and m 2 are both O(l), but note 
that m^, = 0(1/N C ) in the chiral limit - the numerically dominant contribution to its mass 
from the anomaly is formally 1/N C suppressed relative to the 0(1) contribution from the ex- 
plicit chiral symmetry breaking quark masses; the condensates (qq) = 0(N C ); the anomaly 
coefficients a" m = 0(N C ); while finally the coefficient in the topological susceptibility is 
A = 0(1). 

Referring to eqs.(2.11) or (2.12), we therefore see that at large -/V c , x(0) — —A = 0(1). 
Moreover, it is clear from looking at planar diagrams that at leading order in 1/N C , x(0) 
in QCD coincides with the topological susceptibility in pure Yang-Mills theory, x(0)|ym- 
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It follows that 

A = X (0)\ym + O(1/N c ) (2.17) 

a result that plays an important role in the Witten-Veneziano formula (see below). 

Now consider the flavour singlet DGMOR relation (2.7) or (2.14). Each term is 0(N C ) 
apart from A, which is 0(1). Naively, we might think that this sub-leading term would be 
small so could be neglected. However, we know that it contributes at the same order as the 
sub-leading term in (/ 0?? ') 2 Tn^, given by the large anomaly-induced 0(1/N C ) contribution 
to the rj mass squared. So the topological susceptibility contribution to this relation is 
crucial, even though it is formally suppressed in 1/N C . However, we may reasonably expect 
that the further 0(1/N C ) corrections to A are genuinely small and that a good numerical 
approximation to eqs.(2.7) or (2.14) is obtained by keeping only the terms up to O(l). This 
means that we may sensibly approximate the parameter A by x(0)| ym in the DGMOR 
relation. 

This is the crucial simplification. A reliable estimate of x{^)\ym at 0(1) is available 
from lattice calculations in pure Yang-Mills theory [25]. This extra dynamical input allows 
us to determine the four decay constants from eqs.(2.6) and (2.14)-(2.16). We can then 
analyse the flavour singlet decay formula (2.5) to determine gGyy and see how important 
this new term actually is numerically. 

We already know that the Agc^ contribution to eq.(2.12) is suppressed by one power 
of 1/N C compared to the other terms in the formula. Moreover, gc-yy is renormalisation 
group invariant (see refs.[2, 20]). Now, in previous work we have developed an intuition 
as to when it is likely to be reliable to assume that 1/N C suppressed terms are actually 
numerically small. The argument is based on the idea that violations of the OZI rule are 
associated with the U(1)a anomaly, so that we can identify OZI-violating quantities by their 
dependence on the anomalous dimension associated with the non-trivial renormalisation of 
J^ 5 due to the anomaly. RG non-invariance can therefore be used as an indicator of which 
quantities we expect to show large OZI violations. In this case, gc^ is RG invariant, so 
we would expect the OZI rule to be good. This means that since it is OZI suppressed 
(essentially, higher order in 1/N C ) relative to the r( decay coupling <7^ 77 , it should be 
numerically smaller as well. 2 

We will test this conjecture with the experimental data in section 4. The issue is 
an important one. If gcyy can be neglected, then the naive current algebra formulae 
will turn out phenomenologically to be a good approximation to data, even though they 
are theoretically inconsistent. This would apply not just to the radiative pseudoscalar 
decays, but to a whole range of current algebra processes involving the rf (see, for example, 
refs.[19, 20]). We have also used this conjecture in our analysis of the closely related 'proton 
spin' problem in polarised deep inelastic scattering where it plays an important role in our 
prediction of the first moment of the polarised proton structure function g\ [26]. 

Finally, to close this section, we show in detail how these DGMOR relations are related 

to the well-known Witten-Veneziano formula for the mass of the rj' , which is derived in 

2 To be precise, the conjecture is that the contribution ^/GAgc-y-y in the decay formula will be small 
compared to the dominant term f 0r> g n u n . Note that as defined the dimensions of gG-yj and g v '-yy are 
different, so they can not be directly compared. 
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the large- N c limit of QCD. In fact, this is simply the N c — > oo limit of the flavour singlet 
DGMOR formula, which we emphasise is valid for all N c . To see this in detail, recall that 
the Witten-Veneziano formula for non-vanishing quark masses is [6] 

m 2 v , + m 2 v - 2m 2 K = - — X (0)\ YM (2.18) 

Jir 

Of course, only the m 2 , term on the l.h.s. is present in the chiral limit. Now add the 
DGMOR formulae (2.14) and (2.16). We find 

W + W + W + = ^ (2.19) 

To reduce this to its Witten-Veneziano approximation, we impose the large- N c limit to iden- 
tify the full QCD topological charge parameter A with — x(0)|ym according to eq.(2.17). 
We then set the 'mixed' decay constants / 0t? and f Sv ' to zero and all the remaining decay 
constants f 0r i' 7 f Sr i and fx equal to /„-. With these approximations, we recover eq.(2.18). 

In section 4, when we find the explicit experimental values for all these quantities in real 
QCD, we will be able to judge how good an approximation the large- iV c Witten-Veneziano 
formula is to our general DGMOR relation. 

3. Theory 

We now sketch the 1 U(1)a PCAC derivation of the decay formulae and DGMOR relations. 
For a more precise treatment in terms of functional chiral Ward identities and a complete 
renormalisation group analysis, as well as an effective Lagrangian formulation, we refer to 
our original papers [1, 2, 20]. 

The starting point is the U (1)a chiral anomaly equation in pure QCD with n/ flavours 
of massive quarks, viz 3 

d^Jfe = M ab 4 + y^TfQSao (3.1) 

where the axial vector current is J° 5 = q^^b^q and the pseudoscalar is 4% = qjbT a q. 
The corresponding chiral Ward identities for the two-point Green functions of interest are 
therefore (in momentum space): 

i^(J^ 5 Q) - y^daoiQ Q) - M ac (<Pt Q) = (3.2) 

3 We use the following SU(3) notation for the quark masses and condensates: 

(m u \ 

m d = Yl m " Ta 

\ m s ) a=0,3,8 

and 

{uu} \ 

(dd) = 2 

(SS) ) a=0,3,8 

where {<j) a } is the VEV of <ff — qT a q. It is also very convenient to introduce the compressed notation 

M ab = d acb m c $ a6 = d abc {4> c ) 
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ik^(J- 5 4) - V^fSaoiQ 4) - Mac^t 4) = $ab (3.3) 

Since there are no massless pseudo-scalar mesons in the theory, the terms in these 
identities involving the axial current vanish at zero momentum because of the explicit 
factors of k^. The zero- momentum chiral Ward identities simply comprise the remaining 
terms. In particular, we have 

M ac M bd (4 4) = -(M$) ab + (2n f )(Q Q)5 a0 6 b0 (3.4) 

We can also derive the result quoted in eq.(2.11) for the topological susceptibility. In this 
notation, 

X(0) = (Q Q) = — T (3.5) 

1 - {2n f )A(M^ 1 1 

The physical mesons rf = (tt , rf,rf) couple to the pseudoscalar operators 4 and 
Q, so their properties may be deduced from the two-point Green functions above. In 
the PCAC approximation, as explained in section 2, the zero-momentum Ward identities 
are sufficient. In order to make the correspondence between the QCD operators and the 
physical mesons as close as possible, it is convenient to redefine linear combinations such 
that the propagator (two-point Green function) matrix is diagonal and properly normalised. 
We therefore define operators rf and G such that 

(QQ) (Q4)\ ((GG) o 

fiQ (4 4) J ~* V (rfr/) 
This is achieved by taking 

G = Q-(Q4)(( ( f> 5( t> 5 ))ab4 (3.7) 
which reduces at zero momentum to 

G = Q + 2n f A$- b 1 4 (3.8) 

and defining 

V a = f Taa Kb4 ( 3 - 9 ) 
With this choice, the (G G) propagator at zero momentum is simply 

(GG) = -A (3.10) 

and we demand that the (i] a rf) propagator has the canonical normalisation 

(r! a v P) = _ 1 2 S af3 (3.11) 



(3.6) 



k 2 — m 



This normalisation fixes the decay constants introduced in eq.(3.9). The DGMOR rela- 
tions then follow immediately from the zero-momentum chiral Ward identity (3.4) and the 
expression (3.5) for the topological susceptibility. We find 

= -(M^) ab + (2n f )Ad a0 6 b0 (3.12) 
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Unwrapping the condensed SU(3) notation then shows that this matrix equation is simply 
the set of DGMOR relations (2.7)-(2.9). 

The next step is the PCAC calculation of rf — ► 77. We implement PCAC by the 
identification 

<9% a 5 - rml^ + y^- f G5 a0 (3.13) 

which follows from the anomaly eq.(3.1) and the definitions of the fields G and rf . To be 
precise, the — > notation in eq.(3.13) means that the identification can be made for insertions 
of the operators into zero-momentum Green functions and matrix elements only. Notice 
that it is not valid at non-zero momentum, in particular for on-shell matrix elements. It 
is the natural generalisation of the familiar PCAC relation d^J^ 5 — > f^rn^TT, defining the 
phenomenological pion field w. 

The other input is the full axial anomaly for QCD coupled to electromagnetism, viz. 

cKV; 5 = M ab <p b 5 + ^2n^Q5 a0 + a a cm —F^F^ (3.14) 

where a® m are the anomaly coefficients given in section 2. 

Implementing the PCAC relation (3.13) together with the full anomaly equation, we 
therefore find 

^<77l^ 5 |0> = f aa <p(^\0) + 72^7(77|G|0)5 a0 + < m ^( 77 |^^|0) 
= f aa ml^ tf)<mW) + y/*Tf(G G)(77\G}S a0 + < m £( 7 7|^^|0) 

(3.15) 

at zero momentum, where we have used the fact that the propagators are diagonal in the 
r] a , G basis. The l.h.s. vanishes at zero momentum as there is no massless particle coupling 
to the axial current. Then, using the explicit expressions (3.10) and (3.11) for the G and 
r] a propagators, we find the decay constant formulae: 

(X 

/°%«77 + V^/^^^aO = C- ( 3 - 16 ) 

The novel coupling 5^77 is precisely defined from the matrix element (77(G) in analogy 
with eq.(2.3) for the conventional couplings. 

Finally, notice that the mixing of states is conjugate to the mixing of fields. In partic- 
ular, the mixing for the states corresponding to eqs.(3.8) and (3.9) for the fields G and if- 
is 

\G) = \Q) (3.17) 

and 

\v a ) = (r'r ($ab\<t> b 5 ) - ^}A5 aQ \Q)) (3.18) 

In this sense, we see that we can regard the physical rf (and, with SU (3) breaking, the rj) 
as an admixture of quark and gluon components, while the unphysical state \G) is purely 
gluonic. This is why we can usefully picture the unphysical coupling gcyy as the coupling 
of the photons to the anomaly-induced gluonic component of the rj', as already mentioned 
in section 2. 
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4. Phenomenology 



In this section, we use the experimental data on the radiative decays 77, rf — > 77 to deduce 
values for the pseudoscalar meson decay constants f° v , f 0ri , f 8ri and / 8t ' from the set of 
decay formulae (2. 5), (2. 6) and DGMOR relations (2.14)-(2.16). We will also find the value 
of the unphysical coupling parameter gcy-y and test the realisation of the 1/N C expansion 
in real QCD. 

The two-photon decay widths are given by 

3 

rtffoj-ry) = ^lw„l 2 (4.1) 

The current experimental data, quoted in the Particle Data Group tables [27], are 

r(»/-»77) = 4.28 ± 0.19 KeV (4.2) 

which is dominated by the 1998 L3 data [28] on the two-photon formation of the 7/ in 

e + e~ — > e + e _ 7T + 7r _ 7, and 

r(r? -► 77) = 0.510 ± 0.026 KeV (4.3) 

which arises principally from the 1988 Crystal Ball [29] and 1990 ASP [30] results on 
e + e~ — > e + e~7]. (Notice that we follow the note in the 1994 PDG compilation [31] and use 
only the two-photon 77 production data.) 

From this data, we deduce the following results for the couplings gy 77 and <7^ 77 : 

grfyy = 0.031 ± 0.001 GeV" 1 (4.4) 

and 

g vll = 0.025 ± 0.001 GeV" 1 (4.5) 

which may be compared with g^ 77 = 0.024 ± 0.001 GeV. 
We also require the pseudoscalar meson masses: 

Try = 957.78 ± 0.14 MeV m v = 547.30 ± 0.12 MeV 

m K = 493.68 ± 0.02 Mev m w = 139.57 ± 0.00 MeV (4.6) 

and the decay constants f v and fx- These are defined in the standard way, so we take the 
following values (in our normalisations) from the PDG [27]: 

f K = 113.00 ± 1.03 MeV U = 92.42 ± 0.26 MeV (4.7) 

giving f K /U = 1.223 ± 0.012. 

The final input, as explained in section 2, is the lattice calculation of the topological 
susceptibility in pure Yang-Mills theory. The most recent value, obtained in ref.[25], is 

X (0)\ym = — (191 ± 5 MeV) 4 = - (1.33 ± 0.14) x 10" 3 GeV 4 (4.8) 
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This supersedes the original value x(0)|yM — —(180 MeV) 4 obtained some time ago [32]. 
Similar estimates are also obtained using QCD spectral sum rule methods [33]. Using the 
argument explained in section 2, we therefore adopt the value 

A = (1.33 ± 0.14) x 1(T 3 GeV 4 (4.9) 

for the non-perturbative parameter determining the topological susceptibility in full QCD. 

The strategy is now to solve the set of five simultaneous equations (2. 5), (2. 6) and 
(2.14), (2.15), (2.16) for the five remaining unknowns f°v, f 8r i', fv and gGrr . The 

results are 4 : 



that is 



f°V = 10 4.2 ± 4.0 MeV f 0ri = 22.8 ± 5.7 MeV 
f r >' = — 36.1 ± 1.2 MeV f v = 98.4 ±1.4 MeV (4.10) 



J — = 1.13 ±0.04 J — = 0.25 ±0.06 

Jit Jit 

J — = -0.39 ±0.01 J — = 1.07 ±0.02 (4.11) 

Jit Jit 



and 



g Gll = - 0.001 ± 0.072 GeV -4 (4.12) 

These are the main results of this paper. 

Before going on to consider their significance, we can re-express the results for the 
decay constants in terms of one of the two-angle parametrisations used in the literature. 
We reiterate that in our view there is no particular virtue in parametrising in this way, but 
in order to help comparison with other analyses it may be helpful to present our results in 
this form as well. Refs.[14, 15] define 

ff 0r >' /°"\ _ //ocosfl -/osin0 o \ M1 „. 
V/ 8 "' fV ~ \fssm6s / 8 cos0 8 J {4A6) 

Translating from eq.(4.11), we find 

fo = 106.6 ± 4.2 MeV / 8 = 104.8 ± 1.3 MeV 

6»o = - 12.3 ± 3.0 deg 6 8 = - 20.1 ± 0.7 deg (4.14) 

that is 

= 1.15 ±0.05 ^ = 1.13 ±0.02 (4.15) 

Jit Jit 



4 Note that this analysis includes the errors from the experimental inputs and the lattice evaluation of 
x{0)\ym, but not the systematic effect of the approximation A = — x(0)|yM- The errors on the singlet 
decay constants are dominated by the error on A. Isolating this, we have 

f 0v ' = 104.2 ± 0.3 ± 4.0 MeV f v = 22.8 ± 3.5 ± 4.5 MeV 

The octet decay constants are of course unaffected by the value of A. 
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We emphasise, however, that since our definitions of the decay constants differ from those 
in refs.[14, 15], any comparison of the numbers above should be made with care. 

The most striking feature of the results (4.10) for the decay constants is how close the 
diagonal ones, f° v ' and f 8r> , are to /„-, even the singlet. Predictably, the off-diagonal ones 
are strongly suppressed, especially f 0ri . This supports the approximations used in section 
2 in deriving the large N c Witten-Veneziano limit of the DGMOR relations. 

To see this in more detail, it is interesting to compare numerically the magnitudes of 
the various terms appearing in the DGMOR relations, together with their formal orders in 
l/N c . We find (all terms in units of l(T 3 GeV 4 ): 

(/°^) 2 mJ [N c - 9.96] + {f^fm 2 [N c ; 0.15] = [N c ; 0.06] + \f K m\ [N c ; 2.07] 

+ 6A [1; 7.98] 

(4.16) 

f^'f^'m' 2 [N c ;-3A5]+f°V^m 2 [N c ; 0.67] = [Nc - .16] - 2 -^f 2 K m 2 K [N C ;-2M] 

(4.17) 

{f r >') 2 m> 2 [N c ; 1.19] + {f^frn 2 [N c ; 2.90] = - \flm\ [N c ; -0.06] + %f K m 2 K [N c ; 4.15] 

(4.18) 



The interesting feature here is the explicit demonstration that the dominant term in 
the flavour singlet DGMOR (Witten-Veneziano) relation is the topological susceptibility 
factor 6^4, even though it is formally suppressed by a power of 1/N C relative to the others. 
It is matched by the subdominant (O(l)) contribution to (f° v ') 2 m 2 , , which arises because 
of the numerically large but 0(1/N C ) anomaly- induced part of the m 2 , which survives in 
the chiral limit. The numerical results therefore confirm the theoretical intuition expressed 
in section 2. 

To emphasise this point further, we can summarise the numerical magnitudes in the 
combined singlet-octet relation (2.19), which reduces to the full Witten-Veneziano formula 
(2.18): 

(/ "') 2 ^ [N c ; 9.96] + (f^fm 2 [N c ; 0.15] 
+ (f^'fm' 2 [N c ; 1.19] + (f^) 2 m 2 [N c ; 2.90] - 2f 2 K m 2 K [N c ; -6.22] = 6A [1; 7.98] 

(4.19) 

The validity of the large N c limit leading to eq.(2.18) is particularly transparent in this 
form. Numerically, the surprising accuracy of the approximate formula (2.18) is seen to be 
in part due to a cancellation between the underestimates of f s ' n ' (taken to be 0) and fx 
(set equal to /„-). 
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Now consider the decay formulae themselves. The numerical magnitudes and 1/N C 
orders of the various contributions in this case are (in units of 1CP 3 ): 

frfgrfry [N c ; 3.23]+/%^ [N c ; 0.57}+V6Ag Gll [1; -0.005±0.23] = a° cm ^ [N c ; 3.79] 

(4.20) 

and 

/ 8 "V77 [iV c ;-1.12]+/ 8 % 77 [N c ; 2.46] = a 8 m ^ [N c ; 1.34] (4.21) 

The interest here is in the realisation of the 1/N C approximation in the flavour singlet 
decay formula. As explained in section 2, the coupling gcyy is renormalisation group in- 
variant 5 and 0(l/N c ) suppressed and our conjecture is that such terms would indeed be 
relatively small. Remarkably, evaluated at the central value of the topological susceptibil- 
ity found in ref.[25], the coupling gc-yy is essentially zero. This is probably a numerical 
coincidence, since we can not think of a dynamical reason why this coupling should vanish 
identically. What is more reasonable is to consider its value across the range of error of 
the topological susceptibility. In this case, we see from eq.(4.20) that the suppression is 
numerically still under 10%, which is closer to that expected for a typical OZI correction 
although still remarkably small. 

This is a very encouraging result. First, it increases our confidence that we are able to 
identify quantities where the OZI, or leading 1/N C , approximation is likely to be numerically 
good. It also shows that <?g 77 gives a contribution to the decay formula which is entirely 
consistent with its picturesque interpretation as the coupling of the photons to the anomaly- 
induced gluonic component of the rf . A posteriori, the fact that its contribution is at most 
10% explains the general success of previous theoretically inconsistent phenomenological 
parametrisations of rf decays in which the naive current algebra formulae omitting the 
gluonic term are used. 

5. Further discussion 

We now summarise our results so far and discuss a number of aspects of their validity 
and applicability to other processes, both in low-energy pseudoscalar meson physics and in 
high-energy processes such as deep inelastic scattering or photoproduction. 

In this paper, we have seen how the uncontroversial radiative decay formula 

f SV 5V77 + f 8r, 9r)TY = ~ (5-1) 

together with the three DGMOR relations 

(f 0r >') 2 m 2 v/ + (f 0r >) 2 ml = ±(f*ml + 2f*m 2 K ) + 6A (5.2) 

ftff*fm» + f°nf»>m* = ^#(/>4 - SWk) (5-3) 

(/ V)2 m « + (/ a J) 2 m 2 = _l U 2 m 2_± f 2 Km 2 K) (5.4) 
5 The topological susceptibility parameter A is also renormalisation group invariant. 
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lead to the following identification of the four pseudoscalar meson decay constants in the 
mixed rj',rj sector: 

(f* f°n\ / 104.2 ±4.0 22.8±5.7\ 

[f»f H = (,-36.1 ±1.2 98.4±l. 4 J MeV (5 ' 5) 

The novel feature is the use of the flavour singlet DGMOR relation. We emphasise 
again that this is, within the framework of PCAC or chiral Lagrangians, an exact result 
in the sense of being entirely independent of the 1/N C expansion. Its large- iV c limit is the 
well-established Witten-Veneziano formula. This determination of the decay constants is 
therefore on firm theoretical ground and provides a sound basis for phenomenology. 6 

The use of the flavour singlet DGMOR formula relies on the input of the non-perturbative 
parameter A which controls the topological susceptibility in QCD. In time, lattice calcu- 
lations should be able to determine this number accurately from simulations in full QCD 
with massive, dynamical quarks. For the moment, we have to rely on the result obtained 
from the pure Yang-Mills calculation, which we have argued using 1/N C ideas should be a 
good approximation. This temporary approximation is the only place where 1/N C enters 
the determination of the decay constants. 

To illustrate the dependence of the decay constants on the topological susceptibility, 
we have plotted the singlet decay constants / 0?? ' and / 0?? against A in Fig. 1. It is clear 
from the DGMOR and radiative decay relations that the octet decay constants f 8ri ' and 
/ 8t? are themselves independent of A. 





Figure 1: The decay constants f° n and f° v as functions of the non-perturbative parameter A = 
(x MeV) 4 which determines the topological susceptibility in QCD. 

The remaining part of our analysis is the flavour singlet decay formula 

/ 0T? '^ 77 + /°% 77 /°% 77 + V6A 9G77 = ^^EL (5. 6 ) 

Here, because the rj is not a NG boson even in the chiral limit, the naive PCAC formula 

acquires an extra term. In our formulation, this is the parameter gc-yy which we have 

6 Note that the PDG tables currently do not quote values for the flavour singlet decay constants because 
of the subtleties in their definition. A good case can therefore be made for adopting the definitions and 
experimental numbers presented here as a simple and theoretically well-motivated parametrisation of the 
data. 
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argued may reasonably, though certainly non-rigorously, be interpreted as the coupling of 
the photons to the anomaly-induced gluonic component of the rj' , i.e. the component which 
removes its NG boson status. The picture is simple and attractive. The experimental value 
of the unphysical coupling, gc-yy = — 0.001±0.072GeV -4 , means that its contribution to the 
decay formula is under 10%. This is consistent with our expectations for a renormalisation 
group invariant, 1/N C (or OZI) suppressed quantity. 

This is illustrated in Fig. 2, where we have shown the size of the contributions of the 
various terms in eq.(5.6), showing explicitly their dependence on the topological suscepti- 
bility. 



icr 3 




Figure 2: This shows the relative sizes of the contributions to the flavour singlet radiative decay 
formula (5.6) expressed as functions of the topological susceptibility parameter A—(x MeV) 4 . The 
dotted (black) line denotes ^^^f 2 -- The dominant contribution comes from the term f° v gy 77 , 
denoted by the long-dashed (green) line, while the short-dashed (blue) line denotes f av g VT y The 
contribution from the gluonic coupling, y/6AgGjy, is shown by the solid (red) line. 

However, while the flavour singlet decay formula is sensible and theoretically consistent, 
it is necessarily non-predictive. To be genuinely useful, we would need to find another 
process in which the same coupling enters. The problem here is that, unlike the decay 
constants which are universal, the coupling gG-yy is process-specific just like <7,/ 77 or g^-y. 
There are of course many other processes to which our methods may be applied such as 
rj'(rj) — ► V'j, where V is a flavour singlet vector meson p,u>,(j), or rj'(rj) — > 7r + 7r~7. The 
required flavour singlet formulae may readily be written down, generalising the naive PCAC 
formulae. However, each will introduce its own gluonic coupling, such as gcv-y- Although 
strict predictivity is lost, our experience with the two-photon decays suggests that these 
extra couplings will give relatively small, at most O(10 — 20%), contributions if like gc^y 
they can be identified as RG invariant and 1/N C suppressed. 7 This observation restores at 
least a reasonable degree of predictivity to the use of PCAC methods in the U(1)a sector. 

Interestingly, these novel gluonic couplings may also arise in high-energy processes. For 

example, the standard analysis of the two-photon deep-inelastic scattering process e + e~ — > 

7 Notice though that this must be used carefully. For example, the coupling qgnn which enters the 
U(1)a Goldberger-Treiman relation alongside the ry'-nucleon coupling g v 'NN is 1/N C suppressed but not 
RG invariant and is expected to display large deviations from its large N c or OZI limit. See section 6. 
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e + e~ X reduces the problem of finding the first moment of the polarised photon structure 
function gj to a non-perturbative evaluation of the off-shell matrix element (t I =-^5 It) • The 
difference with 7/(77) ~~ * 77 is that in the DIS scenario, the photons are off-shell and the 
interest in the first moment sum rule for gj is precisely how it depends on the target photon 
momentum [7, 8, 9]. Nevertheless, the problem may be formulated in terms of form factors 
grfM-yyik 2 ) and gcry(k 2 ) of which the couplings discussed here are simply the k 2 = limit 
[9]. It has also been suggested that the gluonic couplings gG<tn and gcNN could play a 
dominant role in the photoproduction process jN — > <pN [34]. This interpretation is less 
clear, but it is an interesting subject for future work to look at a variety of electro or 
photoproduction experiments in the light of the PCAC methods developed here. 

6. Pseudoscalar meson couplings of the nucleon and the U(1)a Goldberger- 
Treiman relation 

A further particularly interesting application of these ideas is to the pseudoscalar couplings 
of the nucleon. For the pion, the relation between the axial-vector form factor of the nu- 
cleon and the pion-nucleon coupling g n NN is the well-known Goldberger-Treiman relation. 
Here, we are concerned with its generalisation to the flavour-singlet sector, which involves 
the anomaly and gluon topology. This U{\)a Goldberger-Treiman relation was first devel- 
oped in refs.fll, 12, 13]. In this case, the corresponding high-energy process involves the 
measurement of the first moment of the polarised structure function of the nucleon g± in 
deep-inelastic scattering. In the flavour-singlet sector, this is the so-called 'proton spin' 
problem. (For reviews, see e.g. refs[35, 36].) 

The axial- vector form factors are defined from 



where = 7/7^75 «/27n,jv is the covariant spin vector. In the absence of a massless pseu- 
doscalar, only the form factors Ga(0) contribute at zero momentum. Using the 1 U(1)a 
PCAC substitution (3.13) for d^J^ 5 and repeating the steps explained in section 3 (particu- 
larly at eq.(3.15)), we straightforwardly find the following generalisation of the Goldberger- 
Treiman relation: 



with the obvious definition of the gluonic coupling gcNN in analogy to g v a NN- 
For the individual flavour components, this reads (abbreviating G A (0) = G A ): 



(N\J« 5 \N) = 2m N ( y G A (k 2 )s fM + G a P (k 2 )k.sk^j 



(6.1) 



2m N G a A (0) 



(6.2) 



2m N G A 



(6.3) 



2m N G\ 



f 8V 9-q'NN + f 8V griNN 



(6.4) 



2m N G A 



f° v ' g v >NN + f° v 9 v NN + \/6Ag GNN 



(6.5) 
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Eq.(6.5) is the U(1)a Goldberger-Treiman relation. 8 Notice that the flavour-singlet cou- 
pling G\ is not renormalisation group invariant and so depends on the RG scale. This is 
reflected in the RG non-invariance of the gluonic coupling qgnn [13]. 

In the notation that has become standard in the DIS literature, the axial couplings are 

G\ = \ a 3 G\ = -L a 8 G° A = ±= a° (6.6) 
A 2 A 2^/3 A y J 

and have the following interpretation in terms of parton distribution functions: 

a 3 = Au- Ad 

a s = Au + Ad - 2As 

a = Au + Ad + As-^Ag (6.7) 

Experimentally, 

a 3 = 1.267 ±0.004 a 8 = 0.585 ± 0.025 (6.8) 

from low-energy data, while the latest result for a quoted by the COMPASS collaboration 
[37, 38] is 

«°Iq 2 =4GcV^ = 0.2371°;°;* (6.9) 

It is the fact that a° is much less than a 8 , as would be predicted on the basis of the simple 
quark model (the Ellis- Jaffe sum rule [39]), that is known as the 'proton spin' problem. For 
a careful analysis of the distinction between the angular momentum (spin) of the proton 
and the axial coupling a°, see however refs.[40, 41]. 

From the standard Goldberger-Treiman relation (6.3), we immediately find the follow- 
ing result for the (dimensionless) pion-nucleon coupling, 

9nNN = 12.86 ±0.06 (6.10) 

consistent to within ~ 5% with the experimental value 13.65(13.80) ± 0.12 (depending on 
the dataset used) [42]. 

In an ideal world where g v NN and g v >NN were both known, we would now verify the 
octet formula (6.4) then determine the gluonic coupling gcNN from the singlet Goldberger- 
Treiman relation (6.5). However, the experimental situation with the i] and jy'-nucleon cou- 
plings is far less clear. (See refs.[43, 44] for reviews of the relevant experimental literature 
and recent results.) One would hope to determine these couplings from the near threshold 
production of the r/ and rf in nucleon-nucleon collisions, i.e. pp — ► ppr] and pp — > ppr]', 



8 The original form as quoted in ref.[13] applies to the chiral limit and reads, in the notation of [13] but 
allowing for our different normalisation of the singlet, 

2m N G" A = Fg n , NN H J=F 2 m, 2 „,gcNN 

where F is a RG invariant decay constant defined from the two-point Green function of the pseudoscalar 
field 4> 5 . In the chiral limit, where there is no SU(3) mixing, this is reproduced by the definition (3.12) of 
f 0r> . The off-diagonal decay constant f° v vanishes. The final term is reproduced by eq.(6.5) by virtue of 
the flavour-singlet DGMOR relation (2.7) in the chiral limit. 
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measured for example at COSY-II [45]. However, the rj production is dominated by the S\\ 
nucleon resonance iV*(1535) which decays to Nrj, and as a result very little is known about 
QrjNN itself. The detailed production mechanism of the rj is not well understood. However, 
since there is no known baryonic resonance decaying into Nr/' , we may simply assume that 
the reaction pp —* ppij' is driven by the direct coupling supplemented by heavy-meson 
exchange. This allows an upper bound to be placed on g v 'NN an d on this basis ref.[46] 
quotes g v 'NN < 2.5. This is supported by an analysis [47] of very recent data from CLAS 
[48] on the photoproduction reaction 70 —* prf. Describing the cross-section data with a 
model comprising the direct coupling together with t-channel meson exchange and s and 
-u-channel resonances, it is found that equally good fits can be obtained for several values 
of g-q'NN covering the whole region < g v 'NN < 2.5. 

In view of this experimental uncertainty, we shall use the octet and singlet GT relations 
to plot the predictions for g v NN and gcNN as a function of the //-nucleoli coupling in the 
range < g v 'NN < 2.5. The results are shown in Fig. 3. 



g Bm GeV 3 




Figure 3: These figures show the dimensionless rj- nucleon coupling g V NN and the gluonic coupling 
9gnn in units of GeV -3 expressed as functions of the experimentally uncertain ^'-nucleon coupling 
grj'NN, as determined from the flavour octet and singlet Goldberger-Treiman relations (6.4) and 
(6.5). 

Our main interest in the U (1) a Goldberger-Treiman relation lies of course in the gluonic 
coupling gcNN ■ Unlike its counterpart gcyy in the radiative decay formula, gGNN is not a 
renormalisation group invariant coupling. However, like gc^yy it is suppressed at large N c . 
The various terms in eq. (6.5) have the following orders: G A = 0(N C ), f 0ri , / 0r > = 0(y/W c ), 
A = 0(1), g v NN,g v 'NN = 0(^/TQ, g G NN = 0(1). So the final term Ag GNN is 0(1), 
down by a power of 1/N C compared to all the others, which are 0(N C ). 

The intuition we have developed through experience with flavour singlet physics and 
the large- N c expansion is that while we expect 0(1/N C ) suppressed RG invariant quantities 
to be numerically small, in line with expectations from the OZI rule, we do not expect this 
to be necessarily true for RG non- invariant quantities such as gcNN- 9 So unlike gcy-y in 

9 More precisely, we expect the OZI approximation to be unreliable for quantities which have a different 
RG behaviour in QCD itself and in the OZI limit. The complicated RG non-invariance of gGNN in QCD 
is induced by the axial anomaly, since G° A itself is required to scale with the anomalous dimension 7 of 
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the flavour-singlet radiative decay formula, we would not be surprised if gcNN makes a 
sizeable numerical contribution to the U(1)a GT relation. 10 

This is quantified in Fig. 4, where we have plotted the contribution of each term in 
the U(1)a GT relation as a function of g^'NN- We see immediately that the contribution 
from f° v griNN is relatively constant around 0.08, compared with 2itinG a ~ 0.18. This 
means that the variation of f 0v fly tv n over the experimentally allowed range is compensated 
entirely by the variation of y/QAgoNN- For generic values of g v 'NNi there is no sign of a 
significant suppression of the contribution of the gluonic coupling gcNN relative to the 
others. This should be contrasted with the corresponding plot for 5^77 i n the radiative 
decay formula (Fig. 2). 
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Figure 4: This shows the relative sizes of the contributions to the U(\)a Goldberger-Treiman 
relation from the individual terms in eq.(6.5), expressed as functions of the coupling g V 'NN- The 
dotted (black) line denotes 2itinG a . The long-dashed (green) line is f 0TI g n 'NN and the short- 
dashed (blue) line is f 0r, g n NN- The solid (red) line shows the contribution of the novel gluonic 
coupling, \/6AgcNN, where A determines the QCD topological susceptibility. 

To see this in more detail, consider a representative value, g v 'NN = 2.0, which would 
correspond to the direct coupling contributing substantially to the cross sections for pp — > 
ppr( and jp — ► prf '. In this case, g v NN = 3.96 ± 0.16 and gcNN = —14.6 ± 4.3GeV -3 . The 
contributions to the U(1)a GT relation are then (in GeV): 

2m N G A [N c ; 0.18] = f^' g v , NN [N c ; 0.21] + /°%^^; 0.09] + V6Ag GNN [0(l); -0.12] 

(6.11) 

The anomalously small value of G A compared to G A is due to the partial cancellation of the 
sum of the two meson-coupling terms, which together contribute close to the expected OZI 
value (2m,]\[G A = 0.32), by the gluonic coupling gcNN ■ Although this is formally 0(l/N c ) 
suppressed, numerically it gives the dominant contribution to the large OZI violation in 

the flavour-singlet current J° 5 and all the other terms in the U(1)a GT relation are RG invariant. The 
anomaly, and thus the anomalous dimension 7, vanishes in the large-iV c limit leaving qgnn RG invariant 
as N c — > 00. 

10 Of course, since qgnn is defined with dimension GeV -3 whereas g V NN and g v 'NN are dimensionless, 
we cannot make a direct comparison of the couplings themselves. 
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G A . This is in line with our expectations and would provide further evidence that the 
insights developed in our body of work on both low-energy 77 and 77' physics and related 
high-energy phenomena such as the 'proton spin' problem are on the right track. 

However, it may still be that g v >NN is significantly smaller, implying a relatively small 
contribution to the production reactions pp — ► ppi] and jp — > prj' from the direct coupling. 
In particular, there is a range of g v 'NN around 0.7 — 1.3 where the gluonic coupling gcNN 
only contributes to the U(1)a GT relation at the level expected of a typical OZI-suppressed 
quantity, despite its RG non-invariance. In the extreme case where g^'NN = l-0> we have 
g vNN = 3.59 ± 0.15, g GNN = 0.5 ± 3.8GeV~ 3 , and the contributions to the U(l) A GT 
formula become (in GeV): 

2m N G° A [N c ; 0.18] = f r '' g v > NN [N C - 0.10] + f 0l 'g v NN[N c ; 0.08] + V6Ag GNN [0(l); -0.004] 

(6.12) 

This scenario would be similar to the radiative decays, where we found that the correspond- 
ing coupling 5g 77 ~ using the central value of the lattice determination of the topological 
susceptibility and contributes only at O(10%) within the error bounds on A. This would 
then suggest that RG non-invariance is not critical after all and the 0(1/N C ) suppressed 
gluonic couplings are indeed numerically small. It would also leave open the possibility 
that all couplings of type gcxx are close to zero, which in the picturesque interpretation 
discussed earlier would imply that the gluonic component of the rf wave function is small. 
The suppression of G° A relative to G 8 A would then not be due to gluonic, anomaly-induced 
OZI violations but rather to the particular nature of the flavour octet-singlet mixing in the 
77 — 77' sector. In the parton picture of the 'proton spin' problem, this would be a hint that 
the suppression in a is not primarily due to the polarised gluon distribution Ag but also 
involves a strong contribution from the polarised strange quark distribution As. 

Although we would consider this alternative scenario rather surprising and prefer the 
more theoretically motivated interpretation of the flavour singlet sector in which g v > jy n — 2 
and gluon topology plays an important role, ultimately the decision rests with experiment. 
Clearly, a reliable determination of g v 'NN, or equivalently g^NN, would shed considerable 
light on the U(1)a dynamics of QCD. 
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A. Appendix: Comparison with chiral Lagrangians 

It is useful to compare the results presented in this paper with those arising in extensions 
of the chiral Lagrangian formalism to include the rj and low-energy flavour singlet physics. 
The large- N c expansion plays a crucial role in this approach, since it is only in the limit 
N c — > oo that the anomaly disappears, the chiral symmetry is enlarged to U(3)l x U(3)r 
and the rf appears as a light NG degree of freedom to be included in the fundamental 
fields ip a (a = 0, 1 ... 8) of the chiral Lagrangian. Large- N c chiral Lagrangians have been 
developed by a number of authors, notably ref.[18], though here we shall focus on the 
results obtained by Kaiser and Leutwyler [14, 15, 16]. 

The elegance of chiral Lagrangians should not obscure the fact that all the dynami- 
cal approximations we have made in deriving our results, such as pole dominance by NG 
bosons, weak momentum dependence of pole-free quantities like decay constants and cou- 
plings and a judicious use of the 1/N C expansion, are necessarily also made in the chiral 
Lagrangian formalism, where they are built in to the structure of the initial Lagrangian. 
Indeed, our results can also be systematised into an effective Lagrangian (see ref.[l] for 
details). The principal merit of the chiral/effective Lagrangian approach in general is in 
providing a systematic way of calculating higher-order corrections. 

The physical results obtained using the two methods should therefore be equivalent. 
However, a number of our definitions, notably of the flavour singlet decay constants, differ 
from those made by Kaiser and Leutwyler so the comparison is not straightforward. How- 
ever, we would argue that in many ways the formalism developed in this paper provides 
a better starting point for the description of U(1)a phenomenology, especially in the use 
of RG invariant decay constants and our natural generalisation of the Witten-Veneziano 
formula as the flavour singlet DGMOR relation. 

The fundamental fields in the KL chiral Lagrangian are assembled into matrices U = 
exp[iip a T a ] so that, up to mixing, the fields are in one-to-one correspondence with the NG 
bosons in the chiral and large- N c limits. 11 The Lagrangian is a simultaneous expansion in 
three parameters - momentum (p), quark mass (m) and 1/N C . For bookkeeping purposes, 
KL consider these to be related as follows: p 2 = 0(S), m = 0(5), 1/N C = 0(S), and expand 
consistently in the small parameter 5. It will be clear, however, that this is mere book- 
keeping and should be treated with considerable caution. As we have seen, the realisation 
of the 1/N C expansion in the singlet sector is extremely delicate and it cannot simply be 
assumed that quantities, especially RG non-invariant ones, that are l/N c suppressed are 
necessarily numerically small or indeed of 0(p 2 ,m). 

Nevertheless, arranging the allowed terms in the chiral Lagrangian according to their 
order in 5, KL find: 

C = C + d + C 2 + ... (A.l) 

n We use the same normalisation for the generators as in the rest of this paper, so our singlet field 
differs from the ip of refs.[14, 15, 16] by ip° = ^| ip. The normalisation of the singlet currents and decay 
constants is, however, the same. We assume isospin symmetry and represent the quark mass matrix by 
M = diag(m M , rrid, m s ) with m u = nid = m. 
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where 

£ = -F 2 tr{d^d»U) + ^F 2 Btr{M^U + tfM)-^r{ip°)' 2 (A.2) 

and 

d = 2BL 5 ti(^d ll Wd^U(M^U + WM)^j + 4B 2 L 8 tr ((M^ U) 2 + {W M) 2 ^j 

+ |F 2 Ai^^°a M 99 + i^F 2 BA 2 <p°tr(M j fU - U^M) (A.3) 

Most of the physics we are interested in can be derived from Co + L\. However, we also 
encounter some of the terms in the next order, 

C 2 = 2BL A ti(d^d^U)tT(M^U + U^M) 

+ 4B 2 L 6 (tr(M"<U + WM)Y + AB 2 L 7 {ii{M^U - M)^ 2 

- 2i^BL 18 d^°ti(M j f d»U - d^M) - 4i^B 2 L 2 5<p tr((M j fU) 2 + {WM) 2 ^j 

+ ... (A.4) 

Here, M is the quark mass matrix, B sets the scale of the quark condensate, F is the leading- 
order decay constant before SU(3) breaking, and r is a parameter which is identified at 
leading order in 1/N C with the topological susceptibility of pure Yang- Mills. Their 1/N C 
orders are: F 2 = 0(N C ) and B,t = O(l). The coefficients entering at higher order have the 
following dependence: Ai,A2 = 0(1/N C ), L$,Ls = 0(N C ) and L4, Lq, L7, Lis, L25 = O(l). 

KL define their decay constants Fp (P = ir, K, r], 7/) in the conventional way as the 
couplings to the axial- vector currents: 

(0\J^\P) = ik,F a P (A.5) 

As a consequence, the singlet decay constants F®, and F® are not RG invariant but scale 
with the usual anomalous dimension 7 corresponding to the multiplicative renormalisation 
of the singlet current. It follows that the parameters (1 + Ai) and r are also not RG 
invariant, scaling with anomalous dimension 27. 12 This means that while r coincides with 
A (the non-perturbative parameter determining the QCD topological susceptibility) at 
O(l), it differs beyond leading order in 1/N C . 

The SU(3) breaking which distinguishes the decay constants arises first from the terms 
in C\. Beyond this order, in addition to the direct contributions from the new couplings 

12 Explicitly, the parameters renormalise according to [16] 

TR = Z 2 T B 1+Aifl = Z 2 (l+Ai S ) 

where Z is the usual multiplicative renormalisation factor for the axial current, J^ jR = ZJ%b- in the KL 
chiral Lagrangian formalism, the singlet field <p° is itself renormalised. Allowing for a non-zero vacuum 
angle 6, this field renormalisation is 

& = z-y + yf (z- 1 - 1)9 
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in the chiral Lagrangian, there are also contributions from loop diagrams calculated using 
Cq. These give rise to the 'chiral logarithms' fip = ^n'^F' 2 ^ n ~/7T> wm °h KL have calculated 
explicitly. To present their results, we again use the two-angle parametrisation (cf eq.4.13) 
in the octet-singlet sector: 



t r ] > b r) 



Fo cos i?o —Fo sin i?o 
Fg sin i?g Fg cos t?g 



To this order, the KL decay constants are then [15]: 

4B 



1 + j2 (2^ m 9 L 4 + 2mL 5 ) + 0{pip) 



F K 



Fg = F 



1 + (2 m Q L 4 + ( m + m s)L 5 ) + 0(llp) 

AB / 2 \ 

l + ^[2Y j m q L 4 + -{m + 2m s )L 5 J + 0(jjl p ) 



For the singlet, 



Vl + AiFo 



(A.6) 

(A.7) 
(A.8) 
(A.9) 
(A.10) 



where the scale-invariant part is 

AB 



F 



1 + f2 ( 2 Yl m i Li + 3 ( 2m + m ^)(" L 5 + l a))] (A.ll 



where La = (2L5 + 3Lig)/^/l + Ai = 2L5 + 0(1). Notice there are no loop corrections to 
Fq. Finally, the difference in the angles $0 and $g is determined from 



F o F 8 + F o F 8 



F F 8 sin(i? - 8 ) 
= 8vJB(m-m s )(2L 5 + 3Li 8 ) 



(A.12) 



and is proportional to the S77(3) breaking m, m s mass difference. Also recall [17] that 
at lowest order, the pseudoscalar masses are m 2 = 2mB, m 2 K = (m + m s )B, m 2 = 
|(m + 2m s )B and m^, = |(2m + m s )B. 

These decay constants satisfy a set of relations which are closely analogous, but not 
identical, to the DGMOR relations (2. 14), (2. 15) and (2.16). In fact, up to terms involving 
chiral logarithms, the decay constants shown above satisfy (see also ref.[19]) 



r rf r r,' + r r, r r, 



(F 



Q\2 



(A.13) 



F*F* + F*F* = - F°F 8 sin(i9 - #s) = ^(F 2 - F 2 K ) (A.14) 

F S F S + F 8 F S = (F 8 )2 = _1 {F 2_ 4F 2 k) (A15) 

The corresponding DGMOR relations, i.e. including the appropriate pseudoscalar mass 
terms, are broken also by the terms proportional to L7 and Lg (the Lq contributions cancel). 
These are just the expected 0{m 2 ) corrections to the leading-order PCAC relations. 
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Our main interest, as always, is in the flavour singlet sector. Notice that the relation 
(A. 13) is written for the scale- invariant part of the singlet decay constants only, omitting 
the factor involving the OZI- violating coupling Ai. This does not involve the topological 
susceptibility in any way and is not related to the Witten-Veneziano formula. This enters 
the formalism as follows. In the chiral limit and working to 0(6), the relevant part of the 
chiral Lagrangian is just 

C ~ ^(l + AxJVsV - \r(^f (A.16) 
from which it follows immediately that 



(F^fml, = 6t, F°, = y/T+MF (A. 17) 

At leading order in 1/N C , this is the original Witten-Veneziano formula, since Ai = 0(1/N C ) 
and we may interpret r = — x(®)\ym + 0(1/N c ). However, beyond leading order and away 
from the chiral limit, it does not have a straightforward generalisation in terms of the QCD 
topological susceptibility, since r is not identical to A. Fundamentally, this originates from 
the use of an RG non-invariant field ip° in the formulation of the chiral Lagrangian. In 
contrast, the corresponding formula (2.14), i.e. the flavour singlet DGMOR relation written 
with our RG invariant definition of / 0?? ' and involving the parameter A determining the 
topological susceptibility in QCD itself, is the appropriate generalisation of the Witten- 
Veneziano relation and provides a more suitable basis for describing U(1)a phenomenology. 

These formulae allow a numerical determination of the decay constants and mixing 
angles in the octet-singlet sector. The results quoted in ref.[22] are 

F = 1.25/* F 8 = 1.28/. 
9 = - 4 deg 9 8 = - 20.5 deg (A.18) 

These should be compared with our results, eq.(4.14). Since the definitions of the octet 
decay constants are the same, we expect Fg ~ /g and z?s — #8- Indeed, the angles agree 
while Fg is a little higher than /§, by around 10%. This is explicable by the fact that the 
KL fit incorporates the next-to-leading order corrections in the chiral expansion but not 
the radiative decays, whereas we have used the leading-order DGMOR relation together 
with the octet radiative decay formula. The definitions of the flavour singlet decay constant 
and mixing angle are different in the two approaches, so cannot be directly compared. 

Radiative decays are described in the chiral Lagrangian approach by the Wess-Zumino- 
Witten term, which encodes the anomalous Ward identities. The relevant part, constructed 
by Kaiser and Leutwyler [16], is 13 



aN c 

-wzw - -^r 



tr(eV) + -A 3 tr(e 2 )tr<^ + 2BA 4 ti(e 2 Mip) F^F^ (A.19) 



13 In ref.[16], the constants A3 and A4 are called K\ and K2 resp. The A3 notation is used in refs.[14] and 
[19]. The radiative decay formulae (A. 20), (A. 21) are identical to those quoted in ref.[19], eqs.(39) rewritten 
in our notation. 
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where ip = ip a T a and e = diag(|, — |, — |) is the quark charge matrix. A3 is 0(1/N C ) and 
scales in the same way as Ai, while A4 is O(l) and RG invariant. It is then straightforward 
to derive the following formulae (omitting the 0(m) corrections arising from the A4 term 
for simplicity) 

rfgr/rr + F v9vn = (! + A 3)«em^ (A.20) 

8 8 8 ^ /a 

F-q'9rf77 + FridvTY = a em~ (A. 21) 

These are the analogues of our eqs.(2.5) and (2.6). The key point is that in each case, 
the flavour singlet formula has to incorporate OZI breaking by the inclusion of a new, 
process-specific parameter, A3 in the KL formalism and ga-yy i n our approach, which must 
be determined from the experimental data. In each case, this removes the predictivity of the 
formula unless, as we have discussed, we can bring theoretical arguments to bear to argue 
that the OZI-violating terms are small. In the KL case, this would presumably require the 
RG-invariant ratio (1 + A3)/(1 + Ai) to be close to 1. However, once again, we consider that 
the flavour singlet decay formula presented here has the added virtue of explicitly showing 
the link with the topological suceptibility and giving a physical interpretation of the new 
OZI-violating parameter in terms of the gluonic component of the physical rf meson. 
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